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Abstract 

We present a scalar triplet extension of the standard model to unify the origin of inflation with neutrino 
mass, asymmetric dark matter and leptogenesis. In presence of non-minimal couplings to gravity the scalar 
triplet, mixed with the standard model Higgs, plays the role of inflaton in the early Universe, while its decay 
to SM Higgs, lepton and dark matter simultaneously generate an asymmetry in the visible and dark matter 
sectors. On the other hand, in the low energy effective theory the induced vacuum expectation value of 
the triplet gives sub-eV Majorana masses to active neutrinos. We investigate the model parameter space 
leading to successful inflation as well as the observed dark matter to baryon abundance. Assuming the 
standard model like Higgs mass to be at 125-126 GeV, we found that the mass scale of the scalar triplet 
to be < 0(1O 9 ) GeV and its trilinear coupling to doublet Higgs is < 0.09 so that it not only evades the 
possibility of having a metastable vacuum in the standard model, but also lead to a rich phenomenological 
consequences as stated above. Moreover, we found that the scalar triplet inflation strongly constrains the 
quartic couplings, while allowing for a wide range of Yukawa couplings which generate the CP asymmetries 
in the visible and dark matter sectors. 

Keywords: Cosmology of theories beyond the SM, dark matter theory, leptogenesis, inflation, baryon 
asymmetry, particle physics - cosmology connection. 



1. Introduction 

A widely accepted theory of the early Universe supposes that there has been a period of cosmic in- 
flation fl] |2j [3] which not only explains the drawbacks of standard cosmology, but also provides seed for 
the temperature anisotropy in the cosmic microwave background H [3 Hid El. Finding a particle physics 
model for the inflaton is a non-trivial task however. In the standard model (SM) of particle physics, the 
only scalar field is the SU(2) doublet Higgs , whose quartic coupling Ah is not a free parameter once its 
mass is fixed. Hence a model of chaotic inflation is not possible within the framework of SM. However, by 
adding one more coupling between the Higgs and gravity (9j[l0l[Ill, the potential could be made flat 
enough for producing approximately 60 e-folds of inflation. Indeed there is a plateau for value of the field 
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h » M p i/ ^f^H, where M v \ is the reduced Planck mass. The phenomenological inflationary constraints are 
met when AhI^ 2 h matches the amplitude of density perturbations. For instance with a quartic coupling of 
(9(0.1) the non-minimal coupling to gravity is bounded to be (9(10 4 ), and hence inflation takes place at 
the unitarity scale M pl /£ H - 10 14 GeV 021 El [H EE21 • This is the so-called Higgs inflation |[T6l[T7l[T8ll. 
However, the indication of SM like Higgs at 125-126 GeV 021 [201 lead to a metastable vacuum IT2T1I22I 
at around 10 9 GeV, which is much below the unitarity scale. The current uncertainties in the experimental 
measurements although allow one to extend the vacuum instability up to Planck scale, but it can only be 
resolved at future experiments. One of the possibilities to evade this issue is to widen the scalar field content 
of the SM. Extension of Higgs inflation by means of a scalar singlet or the inert doublet have been discussed 
in (23J M M M Wi- 
lt is paramount to restore a thermal bath at the end of inflation to generate visible and dark matter (DM) 
observed today. At present a number of evidences suggests the existence of DM, which constitutes one 
quarter of the total energy budget of the Universe (251 1291 . However, hitherto a definite mechanism that 
gives rise to the observed relic abundance of DM is unknown. Usually it is assumed that the DM particle 
is in thermal equilibrium in the early Universe and freeze-out below its mass scale (301 • However, an 
alternative scenario to the freeze-out mechanism is that the relic abundance of DM can be accounted by an 
asymmetric component rather than by the symmetric one (3T1 [321 [331 1341 [351 1361 [371 1381 1391 1301 I3T1 1321 1331 

|431|351|3I3|371|3jn|3J[5a[^ Since 
none of the particles in the SM can be a candidate of DM, one needs to explore physics beyond SM to have 
a particle physics candidate for DM. Apart from DM, the non-zero neutrino masses as confirmed by the 
oscillation data are required to be explained in a beyond SM framework. Recall that neutrinos are exactly 
massless within SM because of the conservation of lepton number up to all orders in perturbation theory. 

Besides DM and neutrino mass, an explanation for the observed matter-antimatter asymmetry required 
for the big bang nucleosynthesis is still missing within the framework of SM. If the reheating temperature 
is less than electroweak (EW) scale then it is difficult to generate both DM and the observed baryon asym- 
metry [70] . On the other hand, if the reheating temperature is larger than EW scale, several mechanisms are 
available which can give rise to required baryon asymmetry, while leaving a large temperature window for 
creating DM species observed today. In the past years a lot of effort have been made to unify the mechanism 
giving rise to the asymmetry both in the DM and baryonic sectors OT1 1321 [331 1331 [3511361 [371 [351 1391 1301 13T1 
|321 |4J[3J[3J |371 |MI!2[2a An attempt 

to unify DM and baryon asymmetry via leptogenesis route has also been proposed by two of the authors 
in (581 . where SM is extended by introducing a SU(2)l scalar triplet and a fermionic doublet dark matter 
candidate, stable by means of a remnant Z2 flavour symmetry. The triplet is taken to be at high scale such 
that its out-of-equilibrium decay can produce asymmetric DM as well as visible matter through leptogenesis 
mechanism (761 1771 . Moreover, in the low energy effective theory the induced vacuum expectation value 
(vev) of the scalar triplet could give rise sub-eV Majorana masses to the active neutrinos. Thus a triple 
unification of asymmetric DM, baryon asymmetry and neutrino masses in a minimal extension of the SM is 
achieved. 

In this article, we realize primordial inflation in the presence of non-minimal coupling £a to gravity in 
a scalar triplet (A) extension of the SM and study the consequent low energy phenomenology. An early 
attempt of triplet inflation has been discussed in [78] within the framework of chaotic inflation, where the 
quartic coupling of the triplet is supposed to be negligibly small (less than 10~ 13 ) and the dominant term 
in the scalar potential is the triplet mass, around 10 13 GeV. In presence of the non-minimal coupling of 
the scalar triplet to gravity the mass scale of the triplet can be much below than 10 13 GeV without fine 
tuning the quartic coupling. We take the mass scale of triplet to be around 10 8 - 10 9 GeV such that it not 
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only give neutrino masses, dark matter abundance and baryon asymmetry, but also evade the possiblity of 
having a metastable vacuum in the SM L2UI22- In presence of non-minimal couplings £a and to gravity 
the scalar triplet, together with the SM Higgs field, behaves as inflaton. From this multi-field inflationary 
scenario a single field model can be retrieved as we demonstrate below. We show that once the heavy mode 
is settled down at the minimum, the scalar potential is positive definite only if the mass term and the lepton 
number violating term (jjh^HH) are negligible. However, the inflaton can be an admixture of both triplet 
and SM Higgs moduli or a pure state. We demonstrate in detail how these three cases give rise to different 
constraints on the model parameter space. Subsequently, we explain how the decay of scalar triplet [58] can 
generate an asymmetric dark matter and visible matter observed today. 

The article is organized as follows. In section [2] we briefly underline the main features of the model, 
which has been introduced in ll58l and point out new constraints in the parameter space. We then describe the 
inflationary picture in section |3} where we work out the slow -roll predictions for single field inflation after 
having discussions regarding the numerical and analytical estimates of all the terms in the scalar potential. 
The generation of the asymmetries in the dark and visible sectors are discussed in section [4] The ensuing 
section [5] details the renormalization group (RG) equations accounting for the additional field content with 
respect to the SM ones. Our results are presented in section [6] and we conclude in section [7] We recall 
in | Appendix A|the main Boltzmann equations for the production of the asymmetries in both baryonic and 



DM sectors. 



2. Scalar Triplet as the Origin of Inflation and Darko-Lepto-genesis 

We extend SM by introducing a scalar triplet A(3, 2), where the quantum numbers in the parenthesis are 
the charge under the gauge group S U(2)l X U(l)y- Since the hypercharge of A is 2, it can have bilinear 
coupling to the Higgs doublet H. As a result the scalar potential involving A and H can be given as follows: 

Vj(A, H) = M 2 ^A + ^(A 1 " A) 2 - M 2 H H f H + ^-{H f H) 2 + A& H rfHA*A + — U H A^HH + h.c.l , (1) 
2 2 y2 

where the index J stands for the Jordan frame, as will be explained in the next section [3] The 2x2 
representation of the scalar triplet is 

A-( A+/V " A++ \ (2) 

\ A -A+/V2 " ( ) 



In the fermion sector we introduce a vector-like doublet \p = ((Adm> <A-) with hypercharge Y = -I [58]. 
As a result the bilinear couplings of A to the lepton doublets L, \p and H are given as follows: 

- £ D IjAfDrf + M D M + [f H A l HH + f L ALL + f^A^ + h.c] , (3) 

where fu = ixujM^. The covariant derivative D M is defined as 



D^d^ + iJ-g^ + igatW^, (4) 

where t represents the Pauli spin matrices. For the hypercharge coupling we have used the grand unified 

theory (GUT) charge normalisation: 3 (gf UT ) 2 /5 = (g^f- 
From ([T]) and Q we notice that: 



3 



1. The bilinear coupling of A to the Higgs and lepton doublets jointly violate lepton number by two 
units. Moreover, the couplings are complex and hence can accommodate a net CP violation. As 
a result the out-of-equilibrium decay of A to LL and HH in the early Universe can give rise to the 
observed matter-antimatter asymmetry via leptogenesis route 17511771 . 

2. The Lagrangian is invariant under a remnant Z2 symmetry, with if/ being odd while all the other fields 
even. This ensures the stability of i^dm, the neutral component of if/, which can be a candidate of dark 
matter. Hereafter i^dm is the inert fermion doublet DM [58]. Since the bilinear coupling of A to iff iff is 
in general complex, it can accommodate a net CP violation. Therefore, the out-of-equilibrium decay 
of A — > if/if/ in the early Universe can generate an asymmetry in DM sector in a similar way the lepton 
asymmetry is generated via the decay A — > LL and A — > HH. 

In the effective theory the bilinear coupling of A to HH and if/if/ generates a dimension-five operator 
O5 = if/if/HH suppressed by the mass scale of A. This is an equivalent type-II seesaw for Majorana 
mass of DM. Below EW phase transition this operator generates small Majorana mass for t^DM as 
given by 

m = V2/^(A> = fuU^r , (5) 
Ma 

where v = (H) is the vev of the SM Higgs. Since i^dm is a vector-like Dirac fermion, it can be 
expressed as a sum of two Majorana fermions, i.e. i^dm - (<Adm)l + (^dm)r- Therefore, in a flavour 
basis (((^dm)l» GAdm)d)> the mass matrix of DM is given by 

M = ( M n m 2 ) ■ W 

\m/2 M D ) v ' 

Diagonalising the above mass matrix we get two mass eigenstates (^dm)i and (i/t)m)2 with masses 
Mjj + m/2 and Mp - m/2. The mass splitting 5 ~ m between the two states is required to be (9(100) 
keV in order to explain the high precision annual modulation signal at DAMA lP79l l80l HTH l58l while 
the null result at Xenon 100 [82]. This implies a lower bound on to be 

^ = ^>KT 4 , (7) 
V2(A> 

where we have assumed (A) < (9(1) GeV as required by the p parameter of SM. 

3. In the effective low energy theory the bilinear coupling of A to lepton and Higgs doublets also generate 
a dimension-five operator O5 = LLHH, suppressed by the mass scale of A, for neutrino masses. When 
H acquires a vev, this operator then induces sub-eV Majorana masses to active neutrinos given by: 

r -v 2 
M v = V2/ L (A> - hfn— ■ (8) 
Ma 

For Ma » v, we can easily obtain sub-eV masses of active neutrinos for a wide range of values of 
the couplings fi and For example, taking fi and fu to be order unity we need Ma ~ 10 12 GeV to 
get sub-eV neutrino masses. For lighter A one can get neutrino masses in the ball park of oscillation 
data by taking smaller values of /l, yet maintaining vev of A to be less than (9(1) GeV. An advantage 
for smaller values of fi is that we can easily explain the required ratio: 

R = Ml = tk „ O(10- 5 ) . (9) 
Thus for jfy > 10 -4 , we expect f L > 10 -9 . 



4. In the presence of the non-minimal couplings of A and H to gravity, the scalar potential ([T} can give 
rise to inflation in the early Universe |[T6l [T7l ITSl . The scale of inflation at which the power spectrum 
is normalized (see later section) is [V(A, H)/e] 1 ^ 4 10 16 GeV, which is much below the Planck scale. 
At the end of inflation, the Universe becomes radiation dominated, during which the interactions of 
A as given in (|3]> generate asymmetries in visible and DM sectors. 



3. Scalar Triplet - Higgs Inflation 

3.1. Action in the Einstein frame 

The model for the scalar fields has been defined in the previous section. The scope of this section is to 
work out the action for inflation. The physical fields are defined in the Jordan frame denoted by an index J. 
We introduce for both scalar components non-minimal couplings to the Ricci scalar R. Hence the action in 
the Jordan frame is: 



J d 4 x V=S + + ^A + ex.) R - \<D^H\ 2 - \D M A\ 2 - Vj(H, A) 



(10) 



with the reduced Planck mass set to unity, i.e. M 2 = m 2 j/(87r) = 1. 

In the Jordan frame the couplings £ make the gravitational interactions non-standard. It is therefore 
convenient to perform a conformal transformation into the Einstein frame, for which we put no index, to 
retrieve the standard form of the Einstein equations as far as gravity concern, but at the expense of having 
non-standard kinetic terms for the scalar fields. A conformal transformation preserves the causal structure 
of space-time in both frames and is given by a smooth and strictly positive function of the fields: 

Q 2 = 1 + 2^ A |A| 2 + 2{ H \H\ 2 . (11) 

Note that both frame are equivalent for small field values. The metric and the potential transform as: 

&v = n 2 8^. (12) 

VCft«=™*>. (13) 
The doublet and triplet scalar fields are defined in the unitary gauge as following: 

(14) 
(15) 

where 6 and 6 account for the two degrees of freedom of the triplet neutral component, defined as A = 
[Re(A°) + i Im(A )] / yjl. 

Now taking the large field limit g^d 2 + ^nh 2 » 1 and redefining fields as: 

<p= ^|log(l+^ 2 +^/i 2 ) , (16) 
5 

r = r, (17) 
n 
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( 10 1 reads: 



R 1 



1 r 1 + 1 



2(fo+£ A r 2 ) 3 



2„2 



(5^) - — — - ^ 2 (d II (fi)(d M r) 



,(l_^/V6j 



(«V?) 2 - V(r,p,0) 



(18) 



Note that the kinetic part is highly non-trivial for all fields ip, r and 6. However the potential, with the field 
redefinition, takes the form: 



A H /2 + A HA r 2 + i A r 4 /2 
V(r, </?, 9) = — : — — 1 1 



+ 



A^H + Ur 2 ) 2 



2^H + Ur 2 ? 12 



-2 V /V6 



, ^ + M ^ 2 ,-2,/V6 



2^ H +Ur 2 ) 



-Zp/ V6 ] 



3/2 



(19) 



3.2. Scalar potential analysis 

During inflation the mass eigenvalue of r is very large as compared to the Hubble parameter |[25ll . 
Therefore, r is minimized at ro and we find the effective theory for the light inflatons. The action then 
becomes: 

-2 



1 + 



1+r 2 



(d^) 2 - 



1 



2 ^ + £ A r 2 



(20) 



with V^o, 0) = V(r — > ro, 0). Note that the stabilization of r demands important constraints on the 
couplings, which will be discussed in the following section. For a finite value of ro, with 



, ^-H , 2 ^A 4 

A e ff = — + A HA r + — r . 



we can further approximate the kinetic sector as 

-£kin 



I(^) 2 + i(l-^^)(^) 2 , 



(21) 
(22) 

(23) 



where x = ^o/ Vleff- 

For the potential, as can be seen from ( [T9] >, it consists of three contributions - quartic, quadratic and 
the /i-terms. Since the latter two are exponentially suppressed, one may be tempted to drop them from the 
beginning for simplicity. However we must check explicitly if quartic term is really dominant, only after 
then we can make any simplification. First let us compare the quartic term with the quadratic mass term: 



Vm 



M\r 2 e-^ 



V6&ff 
^eff 



Here we first assume the quartic term is dominant, which normalizes the combination /t e ff/^ 2 



eft 



(24) 



10" 9 from 



the amplitude of the power spectrum (see later section). We will justify this assumption a posteriori. Then, 
with the typical value of (p during inflation, say <p ~ 5, we have e~ 2vl ^ ~ 10~ 2 so that the ratio becomes 



It is not difficult to set this ratio negligibly small with large enough £ e g and not too large r$ and for 
M A ~ 10~ 6 (M A ~ 10 12 GeV), this ratio becomes 10 _5 r 2 /^ e ff which can be easily made small, and even 
easier if we let Ma smaller than 10 -6 . For the triplet term with we can proceed similarly, and obtain 

which looks more stringent than VmIVa and there indeed is a tension: with large enough ro and and 
not too large £ e s this ratio may be close to 1 and we should not neglect V^. However there is another 
constraint that the potential be positive everywhere. For simplicity, let us neglect Vm which can be made 
easily negligible, then the potential is 

V- 10- 10 (l-^^) 2 + ^^cos^^(l-e- 2 ^^) 3/2 , (27) 
which should be positive definite. This gives 

< l0 -^/^^_ e -2 V ,^ l/2 (2g) 

^eff 

1 /2 

We can easily note that eP 1 ^ (l - ^) is a mildly increasing function of <p with the values 1.12364 
at <p — 1 and 7.63495 at tp — 5. Thus, to guarantee the positivity of the potential until the end of inflation 
where tp e ~ 1 provided that V\ is dominant, we should demand 

m^<10- 10 , (29) 

^eff 



l "<10" 13 (^] &. (31) 



which in turn gives, combined with (26 1, 

^~ioV^<kt 2 . (30) 

That is, the positivity of the potential demands that the quartic term be dominant, with the fraction of the 
triplet term contribution at most (9(1) percent. Further, returning back to (25 1, using (29 ) we find 

a 

Vx ~ ~~ \/J.H. 

Thus, for Ma ~ hh, Vm remains indeed negligible compared with V\ unless £ e ff is very large. However too 
large £ e ff will pull down the unitarity scale further, greatly harming the validity of the effective theory: if 
£eff ~ 10 6 > Vm may compete with V\ up to (9(10) percent, but the unitarity scale fiy ~ ~ 10~ 6 may well 
be saturated near Ma ~ Hh ~ 10~ 6 and the low energy approximation cannot be trusted. So not too large 
£ e ft- guarantees negligible contribution of Vm- All these a posteriori justify our assumption at the beginning 
that the potential is dominated by the quartic term so that A e g/^ s ~ 10 . 

This estimate gives us the idea that the contributions of \ to the observable quantities are not significant. 
To check this, we first compute numerically the change in the number of e-folds N as follows. We compute 
N(<p+,x-k) from the moment when the scale of our interest exits the horizon, to e, the end of inflation, with 
a given set of initial conditions of and x* ■ Then we repeats with slightly different initial conditions to 
find N{ip+ + A(p,x+) and Ndp+tX^+Ax). Then, we find 6N according to the change in the initial field values. 
In table [T] we show dN/dtp* and dN/dx* for several values of £ e ff and ro (note that the amplitude of 
the power spectrum fixes Agg for a given £ e ff). Single field analytic estimate ( |52| ) gives dN/dtp* = 35.6967 
(54.0031) for tp* = 5 (5.5). 

7 



Table 1: Numerical estimates of the contributions of the fields <p and x to the total number of e-folds No- For the first case the 
bound \29\ is saturated, while the rest examples satisfy it trivially. 





No 


dN/dy* 


dN/dx* 


r = 1, & ff - 10 4 , n H = 10" 7 , cp* =5, X * = 10- 3 
the same as above but tp* = 5.5 
the same as above but;f* = 10~ 35 


42.0850 
64.3191 
42.0884 


35.9478 
54.2294 
35.9508 


-2.98106 
-5.24260 
-8.89659 


ro = 10, &ff = 10 3 , Hh = 10" y , =5, X * = 10~ 3 
the same as above but (p+ -5.5 
the same as above but^-^ = 10~ 35 


42.1880 
64.5161 
42.1880 


36.0828 
54.4816 
36.0828 


-1.07484 x 10" 3 
-1.98455 x 10" 3 
-3.39071 x 10~ 4 


r - 10 2 , & ff - 50, fi H = 10-", cp* = 5,** - lfr 3 
the same as above but = 5.5 
the same as above but^* = 10~ 35 


42.1785 
64.4986 
42.1785 


36.0711 
54.4600 
36.0711 


-4.19220 x 10- ' 
-9.66338 x 10~ 7 
-3.55271 x 10~ 8 



3.3. Constraints on the scalar potential 

From now on, as discussed in the previous section, we only consider single field case where /j,h term 
does not contribute and only <p drives inflation. However as mentioned at the beginning, we assume that r is 
stabilized already. For this to happen, we need to study in detail this stabilization which gives constraints on 
the couplings. These constraints do affect low energy phenomenology by incorporating RG equations, even 
Hh which does not participate in the inflationary dynamics but whose RG equation does include quartic 
couplings. 

We first have to ensure that the potential, quartic terms alone, is positive definite everywhere. This is 
necessary because we may not have to ensure Aha > 0. Indeed, we must have Ah > 0, Aa > and 

*HA+ V^a >0, (32) 

for positive potential. 

Coming back to the definition of the potential in terms of r and (p, as the mass eigenvalue for r is very 
large compared to H (see Appendix B of [25] and 11831 1841I851 ) we assume r is stabilized at r = ro throughout 
the whole process of our interest. The different minima in which the heavy field r quickly sets in, are found 
minimizing the potential part independent of <p: 

v _ Ah/2 + A A /2r 4 + i tfA r 2 



Vmdep= 4 (fe+ ^ r 2 )2 



The minima are listed below together with the corresponding minimum energy and constraints for vacuum 
stability. At r = and r = 00 inflation is driven by pure Higgs (r = 0) or pure triplet (r = 00). At the finite 
minimum, inflation is driven by an admixture of both fields. 

1. r 2 = (Aha£h ~ Ah£a)/(Aha£a ~ Aa£hY- Then, V^-indep becomes a constant, i.e. vacuum energy, of the 
value: 

VWdep - ^ iXed) - 2 — 7 • (34) 

v p s{A A e H +AHe A -2A H A^H) 
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We demand that v^ mxe ) > and dV 2 /dr 2 \ r2=l 2 > 0. Then, we must satisfy the conditions 



A H A A -A 2 HA >0, (35) 
{hAha-^h <0, (36) 
ZaAha-£hA a <0. (37) 



Note that the first condition is also equivalent to demanding that the numerator of (33 1, which is 
essentially a quadratic equation of r 2 , is always positive, i.e. the equation has no solution of r 2 = 0. 
r 2 — > 0: In this case 8 — > so this corresponds to pure Higgs inflation, i.e. the Higgs moduli alone 
drives inflation. V^indep becomes a constant, i.e. vacuum energy, of the value 



,(H) _ Ah_ 



Vindep = = -| . (38) 



In this case d 2 V/dr 2 \ r 2 =0 > gives 



^ha-^h>0, (39) 
t A A HA - £ H A A <0 . (40) 



r 2 — > oo: In this case h — > so this corresponds to pure triplet inflation (in this case the triplet moduli 
alone drives inflation) with: 

,(A) _ ^A 



Wp^ (41) 



In this case d 2 V/dr 2 \ r 2 =0O > gives 



^hAha-^aAh <0, (42) 
{aAha-Zh*a>0. (43) 



Notice that because of ( 39 1 and ( 43 ) for pure Higgs and triplet inflation Aha > is preferred. 



3.4. Slow-roll analysis for single field inflation 

Provided that the quartic potential alone is dominant over quadratic or triplet contributions to the po- 
tential, we may estimate the inflationary predictions using the so-called SN formalism |[86ll87ll88ll89ll90l . 
Essentially, the 6N formalism tells us that the perturbation in the number of e-folds, which is the same in 
both frames |9"T1 . is equivalent to the curvature perturbation on super-horizon scales. Then the slow-roll 
approximation, described by the parameters e and rj is working welQ 

Before going into the detail of slow-roll inflation let us make comments about the reheating. Inflation 
not only consists of the slow-roll period but also a reheating phase since it permits to link inflation with the 
subsequent radiation dominated era. This phase is connected to the potential part close to the minimum and 
takes place during a few e-folds. The reheating phase is poorly known and technically difficult to model 
properly. To take into account uncertainties on this post inflationary phase we use the reheating parameter 
R m d described in [93 ] as 

AN 

log^rad - +3w reh ) , (44) 



'Note that this approximation is equivalent at first order to the slow-roll predictions obtained with the Hubble flow parameters 
e n , as described in |92|. 
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Figure 1 : Left: Two slow-roll parameters \A1) and {48} are shown as a function of the inflaton field. The magenta solid line denotes 
e while the black dotted line stands for 77, as labeled. Right: The inflationary potential {46} is depicted as a function of the inflaton, 
solid gray line. In both panels the light red region indicates where the slow roll condition e < 1 break down. 



having supposed the simplest model of a scalar field coupled to radiation and that the effective fluid (inflaton 
plus radiation) with energy density p and pressure P is conserved and vv re h stands for the mean equation of 
state parameter w re h = Pip during reheating. In addition AN is defined as the total number of e-folds during 
reheating 

AN = AU - No , (45) 

^V re h being the number of e-folds at which reheating is completed and the radiation dominated period begins 
while A^o is the total number of e-folds during inflation. We assume instantaneous reheating, namely at 
the end of inflation the Universe enters straightaway in the radiation dominated era with equation of state 



w = 1/3. This is equivalent to consider the reheating parameter equal to 1 or log/? ra d = in (44 1. This 
can be understood physically because the pre-/reheating stage can not be distinguished from the radiation 
dominated era and therefore can not affect the inflationary predictions. 

After the analysis of the previous sections the actual potential is, with Vq = A e g /(4^ ff ), 

V(^) = Vo(l-^ 2WVS ) 2 ■ (46) 

Its behavior as a function of the field tp is shown in the right panel of figure [T] Note that inflation takes 
place for trans-Planckian values of the field. The shaded region denotes the breakdown of the slow-roll ap- 
proximation, that we discuss straightaway. We can define the slow-roll parameters in terms of the potential 

as 

1 (V'\ 2 4 g- 4 v/V6 
£ ^[v) =3r—^f' (47) 

77 --— = - V 2 ^V6_L^f . (48) 



V 3 (l_ g -2WV6) 2 



10 



Both parameters, as functions of tp, are shown in the left panel of figure [TJ as labeled. Note that both are 
positive and mono tonic ally increasing functions of the inflaton. The red region indicates the breakdown of 
slow -roll condition e < 1. Then, the number of e-folds becomes, using the slow -roll equation 3H<p + V = 0, 



N 



r* v 

Je V' 



dtp 



V6 



(49) 



To determine the latter, we identify this moment as when e(jp e ) - 1 . Then we easily find 

VS. 



log (2 V3- 3) -0.940. 



(50) 



The total number of e-folds after the Hubble length exit for instantaneous reheating is given by AAf* - 55.6, 
using (44), enough to solve the flatness and horizon problems. We can find <p+ by plugging ip e into (49 1 as 

^ = 5.36. (51) 



From (49 ) we can immediately find 



|^ = ^f e WV6_ 1 \ 4g3 
dtp* 4 V / 



(52) 



The power spectrum is normalized at the pivot scale of WMAP7, Icq = 0.002 Mpc l : Puih)) = (2.43 
0.1 1) x 10~ 9 . The power spectrum of scalar perturbation at the pivot scale is given by: 



V(fP*) ( dN 



Yin 1 [dip*) 2<\n 2 e{ip+) ' 
Thus as quoted several times before, A e g is fixed once we fix £ e ff or vice versa, as 

Meff 



Also, the spectral index is given by 



- 48646.2 V^ff ~ 5 x 10 4 V^ff • 



(53) 



(54) 



fin - 1 - 1e + 1r\ — 



0.965 , 



{dNIdy+Y 

which lies well within the 2cr range of WMAP7, n n = 0.963 ± 0.0014 (29l. 



(55) 



4. Asymmetric DM and Leptogenesis 

From Q we see that there are three different channels available for the decay of scalar triplet A: A — > 
LL, A — > HH and A — » ijjijj. Since the couplings are in general complex, the quasi-equilibrium decay of 
A via these channels produce asymmetries in lepton and DM sectors. The CP asymmetry in either sector 
arises via the interference of tree level with the one-loop self energy diagram as shown in figure [2] From 
these diagrams we see that to generate net CP asymmetry at least two scalar triplets Ai and A2 are required. 
As a result the interaction of Ai and A2 is described by a complex mass matrix instead of a single mass 
term as mentioned in ([TJ. The diagonalisation of the flavour basis spanned by (Ai, A2) gives rise to two 
mass eigenstates = A| 2 Ai + B| 2 ^2 with masses M\ and M2. The complex conjugate of g 12 are given 
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Figure 2: The CP asymmetry for leptogenesis and asymmetric DM generated by the interference of tree and one-loop self energy 
correction diagrams are shown in the first and second raw respectively. 



by £j~ 2 = A~ 2 Ai + B~ 2 A2. Unlike the flavour eigenstates Ai and A2, the mass eigenstates and are 
not CP eigenstates and hence their decay can give rise to CP asymmetry. Assuming a mass hierarchy in the 
mass eigenstates of the triplets, the final asymmetry arises by the decay of lightest triplet ^ and The 
CP asymmetries are defined as 



2 [Br(<f- -» U) - Br(<r+ -» ft)} , 
= 2 [Br (^f -> iAdm<Adm) - Br -> ^dm^dm)] = ^dm , 



(56) 
(57) 



where the front factor 2 takes into account of two similar particles are produced per decay. From figure |2j 
the asymmetries are estimated to be: 



€dm 



1 MiM 2 
8tt 2 M\ - M\ 

1 M X M 2 
8;r 2 M\ - M\ 



where 



Mi 
Mi 

rT 

Mi 



Im 



Im 



Qfi 



2,1, 



(l/ml 2 + l/i/ + I/ilI 2 ) , 



(58) 
(59) 

(60) 



is the total decay rate of the lightest triplet. In the numerical calculations we will use this total decay rate 
as: Ti(m v , Bl, Bh, M\), where Bl and Bh are the branching fractions in the decay channels £1 — > LL and 
fi — > HH respectively. In the following we set m v - 0.05 eV and therefore the total decay rate depends 
only on three variables, namely Bl, Bh and Mi. 

When Ti fails to compete with the Hubble expansion scale of the Universe, fi decays away and produces 
asymmetries in either sectors. As a result the yield factors are given by: 



riL 

Y L = — = e L X ( T] L , 
s 



(61) 
(62) 



where = ng-/s 



Is, s = 2n 2 g*T 3 /45 is the entropy density and tjl, tjdm are the efficiency factors, 



which take into account the depletion of asymmetries due to the number violating processes involving \p, L 
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and H. At a temperature above EW phase transition a part of the lepton asymmetry gets converted to the 
baryon asymmetry via the SU(2)l sphaleron processes. As a result the baryon asymmetry is: 

8« + 4m 



'B = ~ 



I4n + 9m 



Y L = -0.55Y L , 



(63) 



where n is the number of generations and m is the number of scalar doublets. From (62 1 and (63 ) the DM 
to baryon ratio is given by: 

Qdm 1 m DM £dm ?7dm 



Qb 0.55 m p e L T] L 



(64) 



where m p ~ 1 GeV is the proton mass. From this equation it is clear that the criteria Qdm ~ 5 can be 
satisfied by adjusting the ratio 6dm/^l and t]um/^l, where the efficiency factor: 

Yi 



m 



6iX S 



with i = DM, L 



(65) 



can be obtained by solving the relevant Boltzmann equations ||58l |94l |95l given in Appendix A The ratio 
of CP asymmetries is 

£L " Im [(fof^ + flHf* 2H ) Y,aBiflL) a pir 2L )ap] ' 

From the above equation we observe that if > /# » fa, then we get 

epM _ <^(/g) 
e L ~ 0(/ L 2 ) • 

Taking 10~ 5 < fy < 0.1 and /j, ~ 10~ 5 we get the ratio of CP asymmetries in a broad range: 10 2 - 10 8 . 



(66) 



(67) 



5. Renormalisation Group Equations in Scalar Triplet Model 

The RG equations of the scalar, gauge and Yukawa couplings in SM have been extensively discussed in 
the literature, see for example 11961 l97l for discussions relative to the cosmological framework. However, in 
the presence of scalar triplet the RG evolution of these couplings change because of the additional lepton 
number violating interactions of the scalar triplet with SM Higgs and leptons, as it has been described first 
in ||98l and then improved by ||99lll00L which will be our main references. Moreover, in our case, the triplet 
couples to the inert fermion doublet dark matter ijj. In the following we list the modification to the standard 
running for An as well as the RG equations for the different couplings pertaining to A such as A A , Aha, Mh, 

and the non-minimal couplings to gravity £a and 

Having defined fix - dX/d In fi, where fi is the renormalization scale, the RG equations of the quartic 
couplings in the scalar potential including the triplet are given by 

\6n 2 fJ AH =\2A 2 H + 6A 2 HA - ^gj + 9gjj A H + 9 - U^g\ + ^gjgj + g%\ + (\2A H Y? - 121?) , (68) 

16^ = - ^g] + 24gjj A & + + 18^ + ^g\g\ + \4A\ + 4A\ H 

+ 4A A Tr (flf L + flu) - 8Tr (f L f L f L fL + /^/J/,) . (69) 

(9 33 \ 27 
2*1 + y s\\ a ah + + 6g 4 2 + (8A A + 6A H + 4A AH + 61?) A AH 

+ 2Tr (fj L + flu) A AH - 4Tr {f{f L f{fL + f^U) , (70) 
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where we have assumed that the dominant contribution to RG evolution comes from the top quark Yukawa 
coupling Y t in the SM. Note that the beta function for Ah gets a positive contribution from Aha, whose 
importance will be discussed later on. g\, g2 and #3 are the couplings corresponding to the gauge groups 
U(l)y, SU(2)l and S U(3)c of the SM. In presence of the scalar triplet the evolution of gauge couplings are 
given by: 

l^i=^«i. ( 71 ) 

I6n% 2 =-^g 3 2 , (72) 

\6n% 3 =-lg\. (73) 

Since the triplet is a singlet under SU(3)c, the running of #3 is not affected. By the same argument the 
running of the Yukawa coupling for top quark is not affected either: 

l&rfr, = 9 -Y? - {^g\ + 9 -g\ + 8^J Y, . (74) 
The RG evolution of the Majorana Yukawa couplings fi and f^ are given by: 

-3 (flf L + flu) f L -\ \^g\ + 3^J f L + [Tr (fj L + j*f f )] f L , (75) 
I6n%, =3 (flfL + /, - | (Igi + 3glj U + (f L f L + f}fr)] fy ■ (76) 
The RG equation of the coefficient of the trilinear A^HH coupling is given by 

l6n% H = (a h - 4A AH - ^gj - y g 2 + 6F 2 ) m + [Tr (fj L + $fy)] fi H ■ (77) 
The anomalous dimensions of A and H are 

^ 7Ma =|«? + 6g 2 2 - 4A A - Tr(flf L + f*fi) - + . (78) 

2.. _9 2 , 9 o 3, . J, M 2 



16^7M fl + - A 8 Z 2 ~ ~ 37/ - 3 



Mi Mf 1 



where the symbol yx is defined by yx = -X~ l dX/d In fi. The anomalous dimensions have a key role 



in determining the RG equations for the non-minimal couplings and As described in HlOli for a 
general theory of scalars <pi coupled non-minimally to gravity via ^j, the one-loop bare and renormalized 
non-minimal couplings are defined as 

ftfi = (&-^«)^+^y, (80) 

where Z®.. denotes the mass renormalization term m 2 .^. - Z^im^. The RG equations for the £y are linked to 
the mass anomalous dimensions by 

/%; = l6™-^Jyg. (81) 
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Consequently, plugging ( 78 ) into the definition ( 8 1 1, the beta function for the non-minimal couplings are 



16x<P iH = |fe + - \\--gi - -g\ + -A H + 3Yf ) + 3 1 & + - 



2 \ 

Aha + — r 
M 



16/r 2 % =2 |f fl + I J A HA + + I ) 4^ A + I ^ + Tr (/^ + /J/J - ^ 2 _ 6g 



(82) 
(83) 



For renormalization scales below Ma, the triplet is decoupled and should be integrated out. Therefore, 
the set of equations reduces to the SM ones with one important modification. Indeed in the decoupling limit 
(ji < Ma) we see that the trilinear A^HH term provides an effective term 



M 2 



(84) 



As a result the effective quartic coupling of the SM Higgs is modified as 



f t \ 



A - A 



H 



(85) 



with a beta function equivalent to the SM one following the prescriptiorJ^]/l// — > A: 

16^ A = 12A 2 - + A + I ( A g 4 + + + (l2AF 2 _ 12y 4j 

In the decoupling limit, the beta functions of the gauge couplings are also given by 

41 



16;r% =-g\ 



I6n% 2 = 



19 



§2^ 



I6n% 3 =-7gl 



(86) 

(87) 

(88) 
(89) 



The RG equations above and below the mass scale of the triplet are matched at Ma, with particular care 



of the condition in ( 85 1 for the quartic coupling of the Higgs. For the initial conditions at EW scale we 
use the renormalization scale fi = m t as suggested in 111021 . with m t - 172.9 GeV II 1031 and v = 246.22 
GeV. The gauge coupling constants are fixed at the following values: a\{m t ) = 0.01027, a2(m t ) = 0.03344 
and a-i{m t ) = 0.1071. We use the pole matching scheme for /l//(m f ) and Y t (m t ) as detailed in [100] and 
references therein, to relate the physical pole masses to the couplings in the MS renormalization scheme. 
The free parameters Aha, /a» ftp an d Hh are fixed as well at // = m t , the only difference being that their 
running will start only at /u = M A . Below Ma only A will have an effect on the fixed at m t as well. The 
running of the coupling is stopped at the unitarity scale fiy = min(M p i/£//, M p i/£a)- 

The Higgs mass is taken to be 125 GeV, the value at which CMS E0l has reported an excess for a Higgs 
like particle with a significance of about 5cr (and Tevatron HI 0411 with a smaller significance), similarly to 
the 5cr significance of the Higgs boson like found by ATLAS |fT9l at around 126 GeV. The Higgs quartic 



2 We will use both A and A H to indicate the Higgs quartic coupling. It should be intended that in the decoupling limits the two 
are strictly equivalent, while above the triplet mass scale both notations refer to Ah- 
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Figure 3: Evolution of the Higgs quartic coupling A as a function of the renormalization scale ^ until the unitarity scale. In both 
panels, the blue solid line denotes the SM evolution with X H = A and the horizontal gray line denotes the metastability bound. The 
Higgs mass is fixed at M H = 125 GeV and hence A H = 0.26 at the matching scale fi = m,. The common triplet parameters are 
A A = 0.4 and A HA = 0.33, f L = f/, = 0.01 and <f A = £ H = 10 4 . Left, effect of /a h : the green dashed curve stands for \i H = 0.1M A , 
while the magenta dotted line is for \i H = 0.05M A . In both cases the triplet mass is M A = 10 s GeV. Right, effect of M A : the cyan 
dashed line denotes the running for M A = 5 x 10 8 GeV, while the brown dot-dashed line is for M A = 10 10 GeV. Both curves have 
Hh fixed at0.01M A . 



coupling will then have a definite value (Ah = 0.26 at m t ) through the matching pole condition. The running 
of this coupling in the SM is then fixed and goes negative at around 10 9 GeV. A possible way of avoiding 
such a metastable vacuum of the SM is to introduce the scalar triplet at that scale. As a result the positive 



correction proportional to Aha in (68) prevent the Higgs quartic coupling to run towards negative values. 
The role played by the triplet is illustrated in figure [3] The blue solid line stands for the running of Ah 
in the SM: it is clear that at 10 9 GeV this coupling goes below zero, denoted by the gray horizontal line. 
Additionally to this curve, in the left panel we show the running of A for different values of hh taking 



Ma = 10 GeV Since the contribution of //# in (85 I is always negative, the larger its value the sooner 



A becomes negative: the magenta dotted curve is for fy = 0.05, while the green dashed line stands for 
fn = 0.1, which is the upper bound on this parameter in order to prevent instability below 10 8 GeV The 
small step is due to the matching condition at of A and Ah- In the right panel, together with the SM 
running of Ah, we illustrate the role of the triplet mass: the cyan dashed line shows that taking the triplet 
mass at 5 x 10 8 GeV and Aha at least larger than 0.33 will permit the quartic coupling of the Higgs to stay 
positive until unitarity scale, at which it has the value of (9(0.1). On the contrary, if the triplet mass scale is 
too high, as shown by the brown dot-dashed line with Ma = 10 10 GeV, the theory can not be rescued from 
vacuum instability In the above discussions, the values of other parameters can be read from caption. 
Assuming the Higgs mass is fixed at 125 GeV, the following conclusions can be drawn: 

1. Effect of /j,h- The contribution of fin to A is always negative, therefore it lowers the instability thresh- 
old. In order not to loose the interest of the seesaw mechanism to produce neutrino masses in the 
sub-eV range, we keep the mass of the triplet at around 10 8 GeV In addition the ratio fy = fin /Ma 
gives the neutrino mass via type-II seesaw, which can not be negligibly small, hence fy should be at 
least O(l0~ 2 ). We note that those values of fy have negligible effect on the running of Ah and hence 
on the Higgs mass (a large //# could shift the Higgs mass via the tree level relation Mh = VAv) 
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Figure 4: An example of RG evolution as a function fi from m, to the unitary scale fiy. The fixed parameters are: M H = 125 GeV, 
Ah = 0.26, /i = f/, = 0.5, Hh/M& = 0.09. Le/h' Running of the coupling constants in the scalar potential as labeled in the plot. 
Right: RG evolution of the non-minimal couplings to gravity as labeled. 

however may not be too small for neutrino masses. 
2. Effect of In order to avoid instability of the potential, the mass of the triplet should be lower 
than the instability point. Therefore, masses of the triplet larger than 7 x 10 8 GeV are disfavored. In 
the following we fix Ma = 10 8 GeV, which is still in the ballpark to get the correct sub-eV neutrino 
masses, with a corresponding upper bound on fy of (9(0.1). 

In the left panel of figure[4j we show an example of RG running of various quartic couplings in the scalar 
potential, while on the right panel we show the running of non-minimal couplings to gravity, as labeled in 
the plots. Below the mass scale of triplet Ma the quartic couplings Aha and A^ remain constant and start 
to running above Ma due to the globally positive correction from g, and fi. The non-minimal coupling of 
the Ricci scalar with the Higgs receives a negative contribution from the SM parameters below the triplet 
mass scale, while above Ma the corrections due to the triplet parameters set on and increase the value of 
^h, which ends up to be larger than its EW value. On the other hand, £a does not receive any contribution 
from the SM parameters and therefore remains constant up to Ma and then increases. 

We note that for numerical purpose, we use the RG equations with two-loop corrections for SM variables 
while one-loop level corrections for the running of triplet scalar couplings, as described in HlOOi 

6. Results and Discussions 

In this section we discuss the model parameter space which satisfy the inflationary, DM and leptogenesis 
constraints. The following points are in order: 

• The requirement of a positive definite scalar potential leads to a negligible contribution of the angular 
mode 8 for inflation (see table[T]), that is jun < 10~ 7 ; 

• The scalar potential for single field inflation takes the form: V(<p) = vf^ (l - e -2 ^^) , with i = H, 
A and mixed; 
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Table 2: MCMC parameters and priors for the scalar potential parameters and non-minimal coupling to gravity at EW scale fi = m t . 
All priors are uniform over the indicated range. 



MCMC parameter Prior 

Ia 0.1 -» 1 

Aha -1 -» 1.2 

log 10 £• with i = H,A 0^6 

A,/* 0.005 ^1 

ju ff /M A 0.009 -> 0.09 



• The constraints from neutrino masses and Higgs mass (see discussion on and Ma in the previous 
section) sets M A < <9(10 8 ) GeV and fi H /M A = f H <0.\. 

We see that the inflationary constraints are not capable to put stringent limits for the generation of 
the CP asymmetries and vice-versa, because different couplings are involved in each step. We therefore 
discuss separately the inflationary and DM-leptogenesis constraints. The sampling of the parameter space 
is performed via Markov-Chain Monte Carlo (MCMC) techniques, using a modified version of the public 
codes CosmoMC 111051 11061 and SuperBayes B107I . In all plots the triplet mass is fixed at 10 8 GeV and 
M H = \25 GeV. 

6.1. High energy scale — Single field inflation 

The region of the parameter space compatible with inflation is obtained by solving the RG equations 
from EW scale up to unitarity scale (jiu = min(M p i/£#, M p i/£a))> an d by imposing the constraints on both 



the power spectrum measured by WMAP7 (53 1, and the positivity of the scalar potential. The likelihood 
follows a gaussian distribution centered on the measured value of the primordial density perturbations. 
Further constraints are: 



1. All couplings at the unitary scale should satisfy the perturbativity bound: Ai < y4n; 

2. All quartic couplings should be definite positive at all scales, the run down to negative values is 
forbidden: A a, Ah > (below Ma, Ah - A); 

3. Aha + V^ff^A > 0; 

4. The vacuum energy should be positive: V(tp) > 0. 

The sampling is performed over 7 parameters, listed in table|2j together with their uniform prior ranges. 



Mixed inflaton = Vr"l We require additionally the conditions d35t, (136} and (|37l) should be satisfied, 



in order r - r$ to be a positive minimum of the potential V(r), given equation (33). Figure [5] shows the 
ID marginal posterior probability distribution functions (pdf) for all the MCMC parameters at EW scale 
plus the results for Ah at the unitarity scale. The preferred values for A& span all the sampled range with 
a preference for the central values from 0.2 up to 0.8, while Aha is very constrained. This is a result of 
( |35] >. The 2D credible region in the plane MaO^ew), ^#a(j"ew)) is shown in the left panel of figure[6j which 
ultimately indicates that Aha is constrained to be in the range 0.2 — > 0.4 and positive definite. Even though 
the non-minimal couplings to gravity are not observable, we show their preferred value in the third and 
fourth panels of the first raw of figure [5} £a follows a distribution sharply peaked around 10 4 , while £h 
prefers smaller values than £a in a much broader range: from 1 up to 10 4 . Those values prevent the unitarity 
scale to be too low and insure the requirement that the quartic term dominates the potential. The parameter 
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Figure 6: Left: 2D marginal posterior in the j/i A (// EW ), /l# A (p EW ))-plane for mixed inflation. The black solid lines enclose the 
68% and 95% credible region. Center and right: Same as left for pure Higgs and pure triplet inflation respectively. All of other 
parameters in each plane have been marginalized over. 



for lepton violation (second raw, first panel) is loosely constrained: no value is found to be significantly 
preferred. Same for the Yukawa couplings fi and (second and third panels, second raw), except that 
values close to one are slightly disfavored. Eventually the last panel shows Aninu), which is positive 
definite of course and can reach at most values of about 0.3. 



Pure Higgs as inflaton = V Q . We require additionally the conditions (|39| and (|40| should be satisfied for 



the positivity and stability of the vacuum energy. The ID marginalised posterior pdfs are shown in figure [7] 
the first and second panel denote the preferred values for Aa and Aha- These couplings are less constrained 
than the mixed case, since Aha can acquire any positive values from 0.2 up to 1, although the pdf is peaked 
for values of about 0.5-0.6. The 2D credible region as a function of both couplings is plotted in figure [6j 
middle panel: note that it is larger than the mixed case. The non-minimal coupling of the Higgs to gravity 
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Figure 7: Same as figure|5]for pure Higgs inflation. 



is well constrained and its pdf is peaked at 10 4 , while can span a broad range of values from 1 up to 10 4 . 
The parameters related to DM and leptons are unconstrained as in the previous case. The Higgs quartic 
coupling can reach large values up to perturbativity bound. Similar pdfs for all the couplings are obtained 
setting £a to zero. 



Pure triplet as inflaton = Here we require additionally the conditions (|42| and (|43| should be satisfied. 



For all parameters but and the ID marginal posterior pdfs are very similar to the case pure Higgs 
inflation (figure [8]) as well as the 2D credible region for A&, Aha (figure [6j right panel). From figure [8} 
the third and fourth panel (first raw) depict the behavior of the non-minimal couplings to gravity. £h is 
essentially unconstrained and can vary with almost equal probability from 1 up to 10 4 , while £a is described 
by a narrow gaussian centered on its mean value O(l0 4 ), a case similar to mixed inflation. Note that the 
case of £h = gives equivalent results. 

6.2. Low energy scale — DM and visible sectors 



The abundances of matter in the dark and visible sectors satisfying ( 64 1 are produced by the quasi - 
equilibrium decay of the triplet scalar fi. The free parameters are the CP asymmetries e,- for all the species, 
the branching ratios Bj and the dark matter mass hidm- The following constraints apply: 



2> = ' 



\ej\<2Bj. 



(90) 
(91) 
(92) 



The first and third conditions ensure that all amplitudes are physical and the total amount of CP violation 
can not exceed 100% in each channel, while the second one simply demands unitarity of the model. The 
five free parameters are ei, 6dm, Bi, Bdm and widm- In addition to that we allow a hierarchy between the CP 
asymmetries: e# - 6dm ~ 10 2 - 10 5 6l> as remarked in section [4] From RG evolution and vacuum stability 
above 10 8 GeV we require fy < 0.1. Those are the main novelties with respect to the analysis in [58]. 
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Figure 8: Same as figure|5]for pure triplet inflation. 

Table 3: MCMC parameters and priors for the CP asymmetries, branching ratios and m DM . All priors are uniform over the indicated 
range. 



MCMC parameter 


Prior 


log 10 (m DM /GeV) 


1.66 -> 3 


logio ^DM 


-5 -> -2 


togio € l 


-9 -> -4 


l°glO ^DM 


-2^0 


logio B L 


-4.5 -> -3 



The sampling of the parameter space is done via MCMC methods and the parameter inference is 
Bayesian, following the same setup of [58]. We note that the likelihood is given by the sum of a ratio 
function satisfying (64i, that is r = Odm/^> plus a gaussian distribution describing the baryon to photon 
ratio, centered on fjt ± cr^ - (6.15 ± 0.25) x 10"" 10 . The prior ranges are given in table [3] The DM allowed 
range starts at 45 GeV: doublet candidates are excluded below this value by the invisible decay width of the 
Z boson. 

In the left panel of figure [9] we show the ID posterior pdf for rriDM, while all other parameters are 
marginalized over. We see that all the mass range from 45 GeV up to 1 TeV can lead to successful leptoge- 
nesis, namely Yl ~ 10~ 10 and an asymmetric dark matter candidate satisfying ( 64 1. Note from the posterior 
pdf that the most favored region is for low mass candidates, even though there are candidates viable up to 1 
TeV with smaller statistical significance. On the right panel of figure[9]the 68% and 95% credible regions 
are shown in the {twdm> ^DM/^il-plane. From there we see that for DM mass up to around 150 GeV, the 
preferred values of the ratio tjdm/t]l are of O(l0~ 4 - 10~ 5 ), which compensate the large CP asymmetry 
ratio: ebM/^L- However, for DM masses around (9(TeV), eoivi/fL is even larger for the preferred values of 
iJdm/Vl, which decreases down to 10 -6 . For a triplet mass of 10 s GeV the important quantities which drive 
Boltzmann equations are the branching ratios. In figure 10 we show the correlation of t/dm/^l versus Bdm 
and Bi respectively in the left and right panels, within the 68% and 95 % credible regions. We see that the 
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Figure 9: Left: ID posterior pdf for the DM mass otdm- Right: 2D credible regions at 68% and 95% C.L. in the (otdm, '7dm/'7l1- 
plane. All of other parameters in each plane have been marginalized over. 




Io 9io( b dm) Io 9io( b l> 

Figure 10: Left: 2D posterior pdf in the [Bom, ^DM/^il-plane. Right: Same as left in the (Bl, ^dm/'/l) -plane. The credible regions 
are given at 68% and 95% C.L. All of other parameters in each plane have been marginalized over. 



largest efficiency ratio tjdm/^Il is preferred when #dm 1 and small #l — * 10~ 3 - 10~ 4 . This is because of 
the required hierarchy between the sub-eV neutrino mass and the Majorana mass splitting between the DM 
mass eigenstates. We note that without this constraint the preferred values would be the opposite ones, as 
shown in [58]. Since Bum is large, which implies small Bl as 2; Bi - 1 (with i=L, //,DM), the washout is 
large as well, which leads to small ?7dm- On the contrary smaller is the Bl the washout effect is small due to 
inverse decay and hence large ijl is preferred. Note that in either case the production of asymmetry is pro- 
portional to Ti oc 1/ s/BlBh- Therefore when Bl approaches towards 10~ 5 the asymmetry (Yl) as well as 
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Figure 11: Left: 2D posterior pdf in the (?7dm, -plane. Central: 2D posterior pdf in the {B DM , ^ DM )-plane. Right: same as 
left in the [B L ,rj L }-plane. The credible regions are given at 68% and 95% C.L. All of other parameters in each plane have been 
marginalized over. 



the efficiency (t]l) get increased. On the other hand when Bdm approaches towards 1 , which implies small 
Bl, the asymmetry Fdm gets increased but efficiency gets decreased. These behaviors of tjum and tjl can be 
confirmed from figure 1 1 where we have shown the 2D credible regions at 68% and 95% C.L. The extreme 
left one, which constitutes the summary of middle and right ones, reveals that a successful asymmetric dark 
matter and lepton asymmetry can be generated with small t/dm and large tjl- For large B^m and small Bl, 
the CP asymmetry in the DM sector should be larger by more than an order of magnitude with respect to €l 
to compensate the small value of tjdm/Vl- The same behavior is recovered for large DM masses. 

To illustrate the mechanism of the generation of the asymmetries via triplet decay we show two bench- 
mark points in figure 12 which are representative examples from our sampling. In all cases the slow channel 
that builds and conserves the asymmetry is the leptonic one: the smallness of its branching ratio is due to 
the hierarchy between neutrino and DM Majorana masses. The fast channels are both the Higgs, to com- 
pensate the neutrino mass in T\ , and the DM one. Since the DM channel is not related to the neutrino mass 
via T\, its branching ratio can assume different values all along the DM mass range. The first point in the 
parameter space is shown in the left panel, which leads to a successful model for the DM with a mass of 
~ 76 GeV, r ~ 4.7 and Yl - 1.5 x 10~ 10 . The second point in the parameter space is depicted in the right 
panel and accounts for a DM candidate with wdm ~ 60 GeV, r ~ 5.1 and successful baryon asymmetry 
Yl = 1.6 x 10~ 10 . The details about the parameters are given in the caption. For the left panel the branching 
ratios are Bl = 9 x 10~ 5 and Bdm = 0.51, which implies large tjl and small 77dm- Therefore, the ratio of 
Vdm/Vl - 6 x 10~ 6 is small and can be confirmed from figure 10 The fastest channel will be the DM one, 



with the largest branching ratio. For the figure in right panel the branching ratio for leptons is comparable 
with the other benchmark, Bl = 8.5 x 10~ 5 , while the DM one, Bdm = 0.017, is much smaller. This implies 
a larger value for the ratio ?7dm and tjl, because in this case the fastest channel is the Higgs one. The small 
values for the efficiency ratio are compensated by the large CP asymmetry ratio, as already discussed and 
confirmed from figure [lOj 

The DM symmetric component is depleted by the efficient gauge interactions before it freezes out and 
is totally negligible at present day ll56l . while the asymmetric DM abundance in accordance with WMAP is 
proportional to the Yield Fdm- 

Thus we see that in a large portion of the parameter space, in particular around DM masses of (9(100) 
GeV the constraints of having sub-eV neutrino masses and keV mass splitting for mass eigenstates of ifr 
are satisfied. In this case, the ratio of the CP asymmetries ranges from 10 3 up to 10 5 and (67i is easily 
satisfied. We note that those values of the Yukawa couplings are perfectly compatible with the inflationary 
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Figure 12: Absolute value for the Yield of leptons (cyan solid), DM (dotted magenta), Higgs (dashed black), f asymmetry (solid 
red) plus scalar triplet abundancy (black solid) for two successful points as follows. Left: m DM = 76 GeV, B L = 9.1 x 10" 5 , 
Bdm = 0.51, e L = 6.8 x 1(T 8 , e DM = 3.3 x 1(T 4 which leads to r = Sl DM /n b = 4.7, Y L = 1.53 x 1(T 10 and ^dm/^l = 6.49 x 1CT 6 . 
Right: m DM = 60 GeV, B L = 8.6 x 10 5 , S DM = 0.017, e L = 8.9 x lO" 8 , e DM = 1.4 x lO" 5 , Q DM /il b = 5.1 and Y L = 1.64 x lO" 10 
and r] DM /ri L = 2.8 x 10~ 4 . The \Yj\ are rescaled in terms of CP asymmetries. 

constraints. 

Regarding detection constraints for the DM particles, our asymmetric candidate may scatter off nuclei 
in underground terrestrial detectors, giving rise to direct detection signature. Due to the particularity of inert 
fermion DM, the interaction will be inelastic and mediated by the Z boson. For this kind of scattering, a 
Majorana splitting of about 100 keV can explain the DAMA modulated signal and is only partially excluded 
by the upper limit of XenonlOO Il82l released in 2011 and which is the most stringent bound up to now. 
Therefore (Adm as inelastic candidate is allowed in all the range from 45 GeV up to 350 GeV. Further details 
about direct detection of fermion doublet DM can be read from ll58l . 

7. Conclusions 

The indication at the Large Hadron Collider (LHC) of a SM like Higgs boson with mass around 125-126 
GeV suggests that the SM vacuum might be metastable at around 10 9 GeV, although it can be extended up 
Planck scale by considering the present theoretical and experimental uncertainties. In this paper we studied 
the scalar triplet extension of the SM which not only evades the possibility of having a metastable vacuum 
at least up to the unitarity scale, (9(10 14 ) GeV, of the theory but also has a rich phenomenology in presence 
of a vectorial fermionic doublet stabilized by means of a remnant Z2 flavour symmetry and thereby playing 
the role of a DM candidate. 

We introduced non-minimal couplings to gravity for both scalar triplet and the SM Higgs. In presence 
of these couplings the scalar triplet, mixed with the SM Higgs, drives inflation in the early Universe. We 
showed that the extended scalar potential gives rise to slow-roll single field inflation, once the heavy field 
is stabilized at a minimum of the potential. In general the inflaton is an admixture of triplet scalar and the 
SM Higgs. However, depending on the minimum, the inflaton could either be a triplet scalar or be a SM 
Higgs. Taking into account that the potential should be positive definite, these three scenarios give different 
constraints on the quartic couplings, namely A& and Aha- We recall that the quartic coupling Ah is fixed 
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assuming that the Higgs mass is around 125 GeV. 

Unfortunately it is not possible to measure the quartic couplings of the triplet at the LHC, because of 
its large mass. Hence it is not possible to distinguish between the three type of inflationary pictures. Also 
the inflationary scenario does not constrain the Yukawa couplings of scalar triplet to DM and leptons. It 
only constrains the dimensionful coupling between scalar triplet and the SM Higgs to be smaller than 
10~ 7 M p i in order not to destabilize the scalar potential. This is also in agreement with another constraint 
arises from the RG evolution of the quartic coupling of the SM Higgs, which shows that the scale of new 
physics should be order of 10 8 GeV in order not to jeopardize the stability of the SM Higgs potential. Based 
on these constraints we set the mass scale of triplet to be - (9(10 8 ) GeV such that it not only stabilized the 
scalar potential but also gave masses to active neutrinos via type-II seesaw. 

Since the triplet couples to leptons and fermion doublet DM and in general these couplings are complex, 
its out-of-equilibrium decay produce asymmetries simultaneously in either sectors. The lepton asymmetry 
produced by the triplet decay can be converted to observed baryon asymmetry in presence of the EW 
sphalerons. The relic abundance of DM can be accounted by an asymmetric component rather than the 
symmetric component which is usually generated by the freeze-out mechanism. Since DM is a doublet its 
mass is necessary to be larger than 45 GeV in order not to increase the invisible Z width. Moreover, DM 
is inelastic and therefore scattered-off nuclei through Z boson. Since the mass splitting between the two 
companions of DM is about 100 KeV, an order of 100 GeV DM can explain the annual modulation signal 
at DAMA and the null result of XenonlOO. 

An interesting possibility will be to perform a detailed numerical treatment of the inflation scenarios 
taking into account the variation due to the reheating parameter and of multi-field dynamics, on the lines 
of II108II931 . Another attractive possibility is to lower the triplet mass scale down to TeV scale, in that case 
same sign dilepton signal will be accessible at the LHC. The neutrino masses are preserved by means of a 
variant of the type-II seesaw, which involves two scalar triplets II 1091 II 101 . We leave these investigations 
for further works. 
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Appendix A. Boltzmann equations for quasi-equilibrium evolution of triplet scalars 

We briefly report the relevant Boltzmann equations relating the generation of the CP asymmetries in the 
dark and leptonic sectors. A more in depth discussion can be found in [58 ] and references therein. 
The evolution of £f density is described the following Boltzmann equation: 



dX ( 
dz 
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with z - Mi/T and X^ = n^-js = n^Js, if the mass of the triplet stays the same after EW symmetry 
breaking. The decay term is described by 
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while the scattering term and the scattering densities are given by: 
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where H(z) = H(T = M^/z 2 , r = y/l-4/x and x = s/M 2 . 

The asymmetry = (n^- - n^+)/s evolves due to the decay and inverse decay of particles. The 
coiTesponding Boltzmann equation is given by 



dYr 



To 



dz zH(z) 
where Yj = (nj - nj)/s, with j - L,H,\fj and 
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where Xj = nj/s. The evolution of the asymmetries Yj is given by the Boltzmann equation: 



dY 

dz 



k zH(z)X? 



2Y k \ 



(A.11) 



where T^ = y$ /n e ^ is the scattering rate involving the number violating processes, such as LL — > £1 — » 
The front factor in (A. 11 ) takes into account of the two similar particles produced in each decay. Note that 
because of the conservation of hypercharge the Boltzmann equations ( |A.l I, ( |A.9 1 and (A.ll) satisfy the 
relation: 2Y ( + Yj = 0. 
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